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A parametrically forced sine-Gordon equation with a fast periodic mean-zero forcing is considered. 
It is shown that 7r-kinks represent a class of solitary-wave solutions of the equation. This result is 
applied to quasi-one-dimensional ferromagnets with an easy plane anisotropy, in a rapidly oscillating 
magnetic field. In this case the 7r-kink solution we have introduced corresponds to the uniform 
"true" domain wall motion, since the magnetization directions on opposite sides of the wall are 
anti-parallel. In contrast to previous work, no additional anisotropy is required to obtain a true 
domain wall. Numerical simulations showed good qualitative agreement with the theory. 

75.60.Ch 03.40.Kf 



> 

00 

o 

O 

o 

oo 
On 

O 



Oh- 



The sine-Gordon equation (SGE) arises in various 
physical applications including Josephson junction trans- 
mission lines 0J^] , dislocations in crystals || , charge den- 
sity waves || , waves in quasi-one-dimensional ferromag- 
netic materials [DD . The only stable traveling wave so- 
lutions of ordinary SGE for a scalar field <f> are local- 
ized solutions with identical boundary conditions qb — 
and 4> = 27T, called 27r-kinks. On the other hand, non- 
localized kinks separating regions with different values of 
the field (e.g. 7r-kinks) are also important since they are 
easier to follow experimentally. This is especially rele- 
vant for domain wall dynamics in quasi-one-dimensional 
ferromagnets with strong anisotropy of a hard magne- 
tization axis P [7p|JlO[|. In the presence of an external 
magnetic field (weak compared to the anisotropy) ap- 
plied within the easy plane, the spin dynamics can be 
described by SGE (see, e.g., with 27r-kinks as the only 
traveling wave solutions. These localized solutions, pre- 
serving the average magnetization, were observed experi- 
mentally 0Jl^l • An interesting effect of evolution of the 
average magnetization takes place if adjacent domains 
have anti-parallel magnetizations which correspond to n- 
kinks or "true" domain walls. Such 7r-kinks were ob- 
tained by introducing additional anisotropy within the 
easy plane (see, e.g., |J and |jl3), sect. 8.3). 

In this letter we obtain 7r-kinks in a modified SGE with 
a fast mean-zero parametric excitation. For quasi-one- 
dimensional ferromagnets with an easy plane anisotropy, 
this implies the existence of true domain walls generated 
by a rapidly oscillating external magnetic field. As op- 
posed to ||[l3) , no additional anisotropy of the ferromag- 
net is required for the existence of 7r-kinks. 

Various cases of parametrically forced SGE have arisen 



in numerous models of physical systems g,|10|,|lj,|15[ . 

However, the mean-zero case has been left out of consid- 
eration, while it is significant both for theory and appli- 
cations. Indeed, a zero-average periodic forcing implies 
that the states = and (f> = tt are symmetric. This 
leads to the existence of 7r-kinks with such time depen- 
dent macroscopic quantity as the spatial average of the 
field. 

We start with the Hamiltonian P,fl6| 

N 

n = E [ - JS > • S '+! + D (^) 2 - MsHSf ] , (1) 

1=1 

where Sj are dimensionless classical spin vectors, N is 
the number of spins, and J is a constant of exchange 
interaction. The external magnetic field H is directed 
along the X-axis. The anisotropy constant D provides 
the existence of an easy plane, XY, at each site of the 
chain. 

The dynamics of a spin Sj in the effective magnetic 
field is governed by mM 



hSi = -Si x 



dU 



(2) 



Representing spins in spherical coordinates, Si = 
S (cos 9i cos 4>i, cos 9i sin <pi , sin#i), and taking the con- 
tinuum limit, yields a system of two first-order (in time) 
partial differential equations for and <fi |j|l6|. When 
the condition DS ^ g^BH holds ||, the system reduces 
to SGE for the polar angle in the easy plane, cf> 



C 2 



v H sin i 



(3) 
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where C — 2clqS 2 \ // JD / h , v = g/is/^JSaQ , and ao 
is the lattice constant. We take the external field H = 
Hq a(t/e), where a is a mean-zero periodic function with 
a unit amplitude, and Hq is the amplitude of the field. 
Rescaling time and coordinate to dimensionless variables, 
t = tC \JvH^ , x = x \/isHq , we obtain (after dropping 
tildes) a parametrically forced SGE (PSGE) 



4> xx + a(t/e) sin0 = 



(4) 



After averaging Eq. (|]) directly, over the fast time 
scale e, we are left with a linear wave equation with wave 
propagation velocity c = ±1. In what follows, we use 
more subtle averaging to obtain 7r-kink solutions moving 
with any prescribed velocity. 

The phenomenon of 7r-kinks in PSGE has a finite- 
dimensional counterpart: it is the stabilization of the in- 
verted pendulum by periodic vibration of its suspension 
point (Kapitza pendulum lljj]). Since this latter phe- 
nomenon is responsible for the existence of 7r-kinks, we 
outline a very simple geometrical explanation of stability 
of the inverted pendulum with a vibrating suspension. 
Full details with other applications can be found in ]2{J 
along with further references to numerous papers on the 
subject. Consider a pendulum (a bob on a massless rod of 
length V) in a nearly upside-down position, with the sus- 
pension point undergoing vertical vibration. We assume 
the latter to have high acceleration and small amplitude. 
Since the acceleration is large, the force of the rod on the 
bob is large so that the bob would be expected to fol- 
low, in the first approximation, the direction of the rod. 
This suggests considering an auxiliary system where the 
velocity of the bob is actually constrained to the line of 
the rod. In this case the bob will oscillate along an arc 
of a tractrix (the "pursuit" curve: all tangent segments 
from this curve to a straight line have the same length) 
and thus will be a subject to an average centrifugal force 
mk(v 2 ), where k is the curvature of the tractrix and v is 
the speed of the bob which is approximately the speed 
of the suspension point when the pendulum is near the 
top. If we now release the constraint, thus releasing the 
centrifugal force, the bob will behave as if it were subject 
to a centripetal force mk(v 2 ) which pushes the pendulum 
towards the top. If this force exceeds the gravitational 
force, the pendulum is stable; this leads to the simple 
stability criterion (v 2 ) > £g, see [po| . 

To average Eq. (||) we apply a series of canonical near- 
identical transformations via the normal form technique 
jL8|,[l9||, so as to bring the original equation to a better 
form with rapidly oscillating coefficients moved to higher 
order terms. Since the transformations are near-identical 
the solutions for the reduced Hamiltonian are close to 
those for the original one. 

The Hamiltonian of PSGE (g) is given by 



H = 



— a cos < 



dx , 



(5) 



where p = <j) t (below we omit the limits of integration 
and dx) . Let the first canonical transformation be defined 
implicitly as follows 



P = Pi 



SW^p,) 



SW^p,) 
Spi 



The new Hamiltonian is given by H\ = H + W\t , or 



H, = 



a cos ( 



ii 



(6) 
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To kill the rapidly oscillating term in H we choose W\ = 
e Ja_icos0, where a_i is an anti-derivative with zero 
average. The Hamiltonian takes the form 



1 



Pi ^lx ■ • . > 

— + — ea-ip 2 sin </>i + -e a_ x sin q>\. 



(8) 



The last term in the above Hamiltonian cannot be re- 
moved by near-identity transformations since it has non- 
zero mean with respect to t. However, all other terms 
with rapidly oscillating coefficients can be killed. Thus, 
choosing W% — e 2 j a_2P2 sin <\>\ we obtain the Hamilto- 

nian H 2 = J f + % + ±e 2 (a^) sin 2 2 + e 2 R + 0(e 2 ), 
where (R) = 0. Finally, taking W3 = e 3 J i?_i we obtain 
the Hamiltonian 



H 3 = 



Pi 



4, 



1 



sin 



3 + O(e') . (9) 



After rescaling X = ex, T = et, P = 2e 1 p 3 , $ = 2(/> 3 in 
the equations of motion corresponding to (0), we obtain 



Pt 



P + 0{e 2 



xx 



_ x ) sin$ + 0(e) 



(10) 



The system ( |10[ ) is a slightly perturbed SGE with 2ir- 
kinks as approximate solutions. After rescaling back to 
variables ((^3,^3), we obtain 7r-kinks as approximate so- 
lutions 03 w U(x,t), where 



U(x, t) = 2 arctan 



exp e 



x — ct 



(ii) 



Note, that by using the normal form technique, our equa- 
tion (Q) has been brought to the form in which it explic- 
itly represents a slightly perturbed SGE (|Io|). Since we 
have employed only near-identical transformations, the 
original equation (^) should have solutions close to the 
solitary waves given by Eq. (|TT|). 

We have verified our results by the numerical simu- 
lations of PSGE (Q), using the second-order leap-frog 
method. To obtain initial conditions for the original vari- 
ables, 4>(x, 0), p{x, 0) , we start with the initial conditions 
for the transformed variables, 03, p% , generated from Eq. 
(|ll|) as 03(^,0) = U(x,0), p 3 (x,0) = U t (x,0). Retrac- 
ing our canonical transformations and keeping the lowest 
order terms only we arrive at 
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<fi(x,0) = U(x,0) , 

(12) 

p(x,0) = U t (x,0) -ea-i(0)sinU(x,0) . 

In Figure [l] we compare the results ol the simulations 
(dashed line) with the analytical solution ( |TT| ) (thin solid 
line). One can see from the Figure a good quantitative 
agreement between the theory and numerical simulations 
(the two curves in the Figure are almost indistinguish- 
able). We have also simulated PSGE (^), starting with 
initial conditions of shapes different from ( p"2| ) . We have 
found that these solutions split into two linear wave pack- 
ets moving in opposite directions with velocities c = ±1. 
This can be explained by the closeness of PSGE with 
mean-zero excitation to the linear wave equation, as was 
mentioned above. 
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FIG. 1. The behavior of vr-kink solution in the PSGE (|§. 
The thick solid line corresponds to the initial profile given by 
( [i"2| ) with c = 1/2, a = sin(t/e). The thin solid and dashed 
lines correspond to the approximate solution ( |ll| ) and the re- 
sult of simulation, respectively, at t = 300. The analytical and 
numerical curves almost coincide. The parameters of the sim- 
ulations are e — 0.1, time step dt = 0.01, mesh size dx = 0.05, 
system size L = 500. 

In terms of the original physical system described by 
Eq. (|^), the 7r-kink solutions represent moving true do- 
main walls generated by a rapidly oscillating external 
magnetic held. Now we return from the dimensionlcss 
variables to the original variables t = t / C \/vHq , x — 
x I y/isHo , to estimate the values of physical parame- 
ters, which provide the formation of a realistic domain 
wall. The lattice constant clq and energy constants J 
and D are taken, according to Ref. (ll) , as ag ~ 5 A, and 
J ^ D ^ 10 K. In dimensionless units the frequency uj 
of the held and the width of the domain wall Ax are of 



order e 1 . Then, for the held amplitude Ho ~ 1 G and 
e ~ 0.1 0.01 , we obtain lo ~ 10 9 Hz, Ax ~ 10 fim. 
Magnetic held H sin cut, with such values of H and cu 
can be created in an experiment, to observe the predicted 
effect. 

In summary, we have found 7r-kink solutions to PSGE 
with a fast periodic mean-zero forcing. As applied 
to quasi-one-dimensional ferromagnets with the only 
anisotropy of an easy plane, these solutions imply mov- 
ing true domain walls, generated by a rapidly oscil- 
lating magnetic held. Our theoretical results are in a 
good quantitative agreement with numerical simulations 
of the PSGE. These results are also applicable to essen- 
tially two- and three-dimensional easy plane ferromag- 
nets, when plane front solutions are stable with respect 
to small deformations. This problem, as well as the sta- 
bility of 7r-kinks in PSGE, will be addressed in future 
investigations. 

We would like to thank S. Brazovsky, R. Camassa and 
V. Bcilin for fruitful discussions. 
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